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STOCHASTIC PDES: DOMAIN FORMATION IN DYNAMIC 

TRANSITIONS 

GRANT LYTHE* 

Abstract. Spatiotcmporal evolution in the real Ginzburg-Landau equation is studied with space- 
time noise and a slowly increasing critical parameter. Analytical estimates for the characteristic size 
of the domains formed in a slow sweep through the critical point agree with the results of finite 
difference solution of the stochastic PDEs. 

Noise is an indispensable part of the description of symmetry-breaking transi- 
tions, although it is normally included only implicitly in mathematical models. This 
work examines a dynamic transition ( "quench" ) in a spatially extended system. The 
model contains a bifurcation parameter having a critical value at which the state with 
zero order parameter becomes unstable to two equivalent states with non-zero order 
parameter. Above this critical value, the whole system does not choose the same 
state. Rather, a pattern of domains and defects ("kinks") is formed 19, |[ [l6|. Here 



it is shown that the characteristic size of the domains formed when the bifurcation pa- 
rameter is made explicitly time-dependent depends logarithmically on the magnitude 
of additive noise: additive noise is thus a determining influence even when extremely 
small. The effect of low-intensity multiplicative noise is much less dramatic. Noise is 
added in such a way that it has no correlation length of its own (white in space and 
time). Numerical results are reported, obtained using a finite difference algorithm 
whose continuum limit is an SPDE. The simplest model with the essential features, 
the real Ginzburg-Landau equation, is considered here. The dynamics take the system 
from a potential with a single minimum to one with two minima (Figure 1). Similar 
results have been found for the Swift-Hohenberg equation, which is more explicitly 
designed to model Rayleigh-Benard convection jl9) . 

The Ginzburg-Landau SPDE. The order parameter Y t (x) is a function of space 
and time and the bifurcation parameter g. For g < fixed, the state {~Y t (x) = Vx} 
is stable, and only noise prevents the system from approaching it arbitrarily closely. 
For g > fixed, one sees a pattern of regions in which Y t (x) is positive and regions 
in which Yj(x) is negative (domains) separated by narrow transition layers (kinks). 
The subject of this paper is not the slow merging of domains or the nucleation of new 
kinks - these happen on timescales much longer than those considered here [|. 
This paper derives the number of domains formed when the parameter g is explicitly 
time-dependent, starting from g < and ending with g > 0. 

The SPDE describing the dynamics is written in the following dimcnsionless form: 

(0.1) dY t (x) - {g(t)Y t (x)- Y t (x) 3 + CY f (x))dt + edW t (x). 



l l 
i' \ i' 

Brownian sheet. The equations were solved as initial value problems, with 



Here Y t (x) : [0,L] m x [— i, i] x £1 — > 1Z, ft is a probability space and W t (x) is the 



(0.2) g(t) = /it 



* Centre for Nonlinear Studies, Los Alamos National Laboratory, New Mexico NM87545, USA 

1 



2 



Grant Lythe 



slowly increased from —1 to 1. Periodic boundaries in x are used so that any spatial 
structure is not a boundary effect. The constants /i, e and \ are all <C 1. The spatial 
operator C — A where A = Y^iLx ~§~7i the Laplacian in lZ m . 

An alternative scaling of ( |0.l[ ) is sometimes illuminating: if x is rescaled so that 
[0, L] -> [0, 1] then Q becomes 

(0.3) dY t (x) = {g(t)Y t (x) - Y f (x) 3 + DAY t (x)) dt + D^edW t {x), 
where D = L 2 . 

The Brownian Sheet. To add noise to an ordinary differential equation, one adds 
increments of the Wiener process at each time step; the idea is of a particle contin- 
uously subject to small impulses ||. To add noise to a partial differential equation 
(PDE) one adds increments of the Brownian sheet, that are random in space and 
time |24|, |2(|. One can imagine, for example, the motion of a flexible sheet in a sand- 
storm. Typically the motivation for adding noise to PDEs is to make allowances for 
small-scale, rapidly-varying effects with mean zero that have been neglected in the 
derivation of an equation. 

The Wiener process can be thought of as assigning to each successive interval of 
the time axis a Gaussian random variable with variance proportional to the length 
of the interval. The Brownian sheet assigns to each volume element in 1Z + x lZ m a 
Gaussian random variable whose variance is proportional to the volume of the element 
p4f . More precisely, it is possible to define a map A from B(H+ x lZ m ) to a probability 
space such that for each h E B(TZ+ x 7Z m ), A(h) is a Gaussian random variable with 
mean zero and (A(hi)A(h 2 )) — l{h\ n h 2 ) where I is the Lebesgue measure Q. The 
Brownian sheet, so called because its realisations in m — 1 look like a ruffled bed-sheet 
tucked in on two adjacent sides, is defined as W t (x) = A([0,t] x [0, #]), where [0, ir] 
is the element (interval, square, cube, . . .) with opposite corners at the origin and at 

x e n m . Thus 

(0.4) W t (x) : n x K+ x TZ m -> TZ. 

The set is the set of labels for realisations; averages over realisations are denoted by 
angled brackets. Each W t (x) is a real- valued Gaussian random variable with mean 
zero and variance (Wf (x)) — tx m . Stochastic processes will be denoted in bold with 
the subscript t. For example, the Wiener process is denoted Wf and is constructed 
by defining W t = A([0,t}). We denote processes satisfying SPDEs with upper case 
bold Roman letters with subscript t and x in brackets. 

The action of Wt(x) as an integrator is easily described in the case where the 
integrand is a deterministic function. If fx(x,t) and f2(x,t) are continuous functions 
on V x [0, t], where V S K"\ then 

It= / [ fi(x,s)dxdW s (x) and J t =[ f f 2 (x, s) dx dW s (x) 
Jo Jv Jo Jv 

are Gaussian random variables with < I t >= 0, < Jt >= and 



< It J t >= / / f 1 (x,s)f 2 (x,s)dxds. 
Jo jt> 
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Numerical Solution. An approximate numerical solution of an SPDE is produced 
by solving a large set of SDEs. The finite difference method for a parabolic SPDE 



consists of replacing the infinite dimensional system (3.1) by N m ordinary SDEs on 
a grid of equally-spaced points in [0,i] m separated by Ax. The SDE at site i is 

(0.5) dY t (i) = {g(t)Y t (i) + YUi)) dt + (Ax)" 2 AY t (i) dt + (Ax)-^ 2 edW t {t), 

where the discrete Laplacian A is defined by 

(0.6) AY(i) = Y ( i ') - 2mY(i) 

i' 

and the sum is over the 2m nearest neighbours of x. The random variables added 
at neighbo uring grid points are independent no matter how small Ax is. Stochastic 
PDEs like fl0.l|) use space-time white noise along with a deterministic operator that 
includes some spatial coupling (here the Laplacian) to generate fields with non-delta- 
function correlations in space and time. Linear SPDEs can be very efficiently solved 
in Fourier space |t], However, the finite difference method is more easily adapted 
to nonlinear SPDEs. One finite-difference numerical realisation of an SPDE generates 
an approximate solution for one lo € f2, at a discrete set of times {ti} and positions 
{xi}. The time-stepping used to generate the figures in this work is the stochastic 
analogue of the second-order Runge-Kutta method jl0| [l4|, |l3) . 

Fourier decomposition. For k 6 Z rn , let the Fourier coefficient yt(fc) be defined 

by 

(0.7) y t (k) = L~^ ( e iKX Y t (x)d m x where k = — k. 

J[o,l]™ L 

Each yt(k) is a complex- valued stochastic process satisfying the following SDE when 
k > 1: 

(0.8) dy t (k) = (g(t) - n 2 )y t {k) dt + dW t (fc) + nonlinear terms. 

v2 

Since Y t (x) is real- valued, for k < — 1 we obtain yt(fc) from the relation yt(— k) = 
(y t (k))* . The imaginary part of y(0) is always 0; the real part satisfies 

(0.9) dy t (0,r) = g(t)y t (0,r) dt + edW t (0,r). 

Each Wt(fc) is an independent complex- valued Wiener process. When |y*(fc)| <C 1 
for all k, the nonlinear terms that involve products of yt(fc)s are unimportant and the 
SDEs ( |0.8[ ) can be solved separately. It then becomes necessary to evaluate the value 
of g at which these nonlinear terms become important. This is called an exit value 
problem in a dynamic bifurcation |l^, |2l], |2^, [2j| . 

Dynamic bifurcation. We consider a stochastic ordinary differential equation slightly 
more general than ( |0.9| ), including also multiplicative noise: 

(0.10) dy 4 = g(t)y t dt + edw t { V + 7 y t dw t ( 2 \ 

where w^ 1 ' and are independent real-valued Wiener processes, g = fit and the 

initial condition is y t = yo at g = pio = — 1. < e, 7 <C ^f\l -C 1. The solution of 



is 



y, = !ln rxp ( ^(f 2 - t 2 ) - i 7 2 < + 7W t < 



(O.llf e exp Q^ii 2 - ^ 7 2 t + 7w ( (2) j exp ^-^A*s 2 + ]p 2 s - 7 w s 



( 2 M dw w 
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Fig. 0.1. Dynamic transition in terms of a potential. The dynamics take the system from 
a single-well potential when the critical parameter g is negative to a double-welled potential for g 
positive. When g is explicitly time- dependent, the system lingers near the origin until well after 
g = 0. When the increase of g is sufficiently slow, the value of g when the system finally moves 
away from the origin is related to the noise level e and the rate of increase of /i by g = yj 2fi log e. 



The second term in ( 3.11 ) has zero mean and variance given by 



(0.12) (y t 2 ) - (y t ) 2 = e 2 exp(/it 2 + j 2 t) / exp(- M s 2 - 7 2 i) ds. 



to 



For g > y/JI, this variance is well approximated by 



(0.13) 



2\ 



y* 



Comparing the mean and standard deviation of yt for g > we find that, for 
fi, 7 <C 1 and \yo\ < C(r)i the standard deviation is much larger than the mean if 
/i \ loge| < i. Then the dynamic bifurcation is noise-dominated Jig] , y t remains small 
until well after g passes through 0, and finally attains an 0(1) magnitude at a value 
of g with mean value g = g c + O(fi) and standard deviation proportional to /i where 



(0.14) 



9c 



^|loge| 



Stochastic equations like ( JO. 1| ) and ( 0.10| ) with g constant and 7 7^ have attracted 
much attention because multiplicative noise acts in a way that can be interpreted 
as a shift of the bifurcation point from g = j2|. However, in all cases the shift 
is proportional to 7 2 compared and is dwarfed by the large dynamic delay of the 
bifurcation of interest here. In (3.13) we see that, for the variance of yt, g is replaced 



by g + • Note, however, that for 7 ^ the distribution of yt(x) is not Gaussian. 
A slight non-Gaussianity is the main effect of multiplicative noise in the SDE and 
SPDE systems of interest here and will not be further discussed. 

Space-time evolution. Now return to the spatially extended system. We first 
make the hypothesis that, as in the dynamic bifurcation, Y t (a;) remains everywhere 
small for g < g c = \/2fi\ loge|. Then the emerging pattern of domains can be studied 
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from the linearised version of (0.1) (i.e. without the cubic term). The solution of the 
linearised version of (|0.l|) is: 



(0.15) Y t (x) = / G(t,t ,x,v)f(v)dv + e f f G(t, s, x, v) dv dW s (v) 

J[0,L] m Jto ^[0,L] m 

where G(t, s, x, v) — 

(0.16) (4 7 r(t- S ))-/2 ex p(-^_ s 2 ) ) J2 eip M'^f 



j=-oo 



A(t - s) 



The first term in (0.15), dependent on the initial data /(x), relaxes quickly to very 
small values. The correlation function is therefore obtained from the second, stochas- 
tic, integral in (0.15). Performing the integration over space, assuming that L > (-)2, 
gives 

(0.17) c(x) = (Y t (x')Y t (x' + x)) = e 2 / — — _- ds. 



la 



in(t-s)Y 



Sincc Yj(x) satisfies a non- autonomous SPDE, the correlation function is explic- 
itly a function of time. In the early part of the evolution, however, the deviation 
from that obtained from the corresponding static (<? =constant) equation is small. 
We consider this quasi-static period by making the change of variables u — t — s in 
Then 

(0.18) exp(/x(< 2 - s 2 )) = exp(2fitu - fiu 2 ) = exp(2^iu)(l - fiu 2 + ...). 
and 

£ 2 / roc e -2\g\u _^ roo e -2\g\u _^ 

(0.19) c(x) = —— — ^ / j^-e 8"da-/i jrr-^e ^ du + . 

^ ' W (8tt)" 1 /2 yj Q u m/2 r J q u m/2-2 



e 2 \g\ * 2 i_m / r , / /t-Tn A 6 ^ 2 



Thus 

(0.20) c(x) = -ig-^x-^ ^^(xVlsl) + ^\g\K^ 3 (xV\g\) 

where K m is the modified Bessel function of order m. For example, when m = 1, 
(0.21) c(x - s') = -^e-l-^lVM ( 1 - 14(1 + x^\g~\ + \x 2 \g\) + 



The first term in (0.20) is the (long-time) correlation function for the SPDE obtained 
by fixing g < JL5| . Clearly the expansion in p- is no longer useful for g > —y/Ji. 
The correlation function itself, however, remains well-behaved as g passes through 0; 
the divergences associated with critical slowing down are not present. 

In one space dimension, the solution of the SPDE ( p.l[ ) is a stochastic process 
with values in a space of continuous functions [ p4| That is, for fixed u> £ Q and 
t G [— i, !■], one obtains a configuration, Yt(x), that is a continuous function of x. 
This can be pictured as the shape of a string at time t that is constantly subject 
to small random impulses all along its length. In more than one space dimension, 
however, the Yj(x) are not continuous functions but only distributions |24| , ^| and 
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Fig. 0.2. Mean squared value of Yt(x). For g > ^/jS, tfte variance o/Yt(x) is prop ortion al to 
cxp(/it 2 ). T/iis exponential rise eventually takes ~Yt(x) to 0(1) values. The solid line is (0.2i) and 
the dots are from a numerical realisation in two space dimensions with fi = 0.002, e = 10~ i0 and 
L = 12tt, carried out on a 256 X 256 grid. The vertical dotted line is at g = g c = \J 2p\ log e| . 



the correlation function c(x) diverges at x = 0. In the non-autonomous equations 
studied here, however, the divergent part does not grow exponentially for g > 0, and 
by g = y/2fi\ log e\ it is only apparent on extremely small scales, beyond the resolution 
of any feasible finite difference algorithm. 

We now examine the evolution for g 
using Laplace's formula: 



> 



yZ/z, where we can approximate (0.17) 



(0.22) 



c(x) 



3 fit 



V/i (8i)W2 



Thus typical values of Y t (x) increase exponentially fast and the correlation length at 
time t is y/St. Once c(0) > 0(e), the noise no longer greatly influences the evolution; 
its effect can be thought of as wiping out the memory of the initial condition at g < 



and replacing it with an effective random initial condition. From (3.22), we see that, 
at g = y2/x] loge|, c(Q) = 0(1) and the cubic nonlinearity can no longer be ignored. 
(Figure 2.) 

Density of kinks. The second hypothesis that enables the characteristic domain 
size to be estimated is that the effect of the cubic nonlinearity, when it finally makes 
itself felt, is to freeze in the spatial structure. This is indeed the case in numerical 
simulations: no perceptible changes occur between g = g c and g = 1 (l7[ [l9|]. Thus 
the correlation length at g — g c , 



(0.23) 



A = = 2 



_ 9 7/4 



iogfl 



becomes the characteristic length for spatial structure after g = g c . In one space 
dimension it is possible to put the scenario just described to quantitative test by 



SPDEs 



producing numerous realisations of the non-autonomous SPDE ( |0.l| ) and recording r, 
the number of times that Y t (x) crosses upwards through in the domain [0, L] at 
9 = 1. The average can then be compared with the mean number of upcrossings for 



Gaussian random field with correlation function (0.17) at g = g c . 

Consider a homogeneous Gaussian random field Y t (x) with correlation function 
c(x). Then Y t (x + Ax) — Y t {x) is a Gaussian random variable with mean zero and 
variance ((Yt(x + Ax) — Yt(x)) ) = b(Ax), where 

(0.24) b(Ax) = 2 (c(0) - c(Ax)) = 2c'(0)Ax + c"(0)Aa; 2 + . . . . 

The probability that Y t has an upcrossing of between x and x + Ax is given by 

^[(upcrossing S (x,x + Ax))] 

[ V[Y t {x) = u}V[Y t (x + Ax) - Y t (x) > u]du 

J —OO 

=(2vr)- 1 (&(Ax)c(0))-^ / e-" 2 / 2c (°) / c^W^&v du 



o 



1 fb{Ax)Y f°° . 2 b(Ax).._ i ,, , 
1 1 ' exp(-w^ v ; )(l-erf(w))dw 



2 V ttc(O) J J ^ v c(0) 

1 / fb(Ax) 

(0.25) = — arctan ' 1 



2tt U c(0) 

Now consider a grid of total length L made up of N sites separated by Ax. Let 
Ax — > with L fixed, i.e. let TV — > oo. When c'(0) 7^ 0, the number of zero crossings 
of Y t defined on this grid is proportional to As"! for Ax — > 0. When, as in ( p,17| ), 
c'(0) = 0, the mean number of zero crossings per unit length approaches a finite 
number as Ax — ► given by El Q : 




(0.26) 

In Fig. 3 this average is displayed as a function of the sweep rate for the case ( |3.22] ), 
evaluated at g = g c : 

(0 - 27) r ^ h 



Summary. When the critical parameter in the Ginzburg-Landau equation is slowly 
increased through a characteristic length is produced as follows. The field remains 
everywhere small until well after the critical value for the loss of stability of the uni- 
form symmetric state. Meanwhile, the correlation length grows proportional to \ft. 
At g ~ g c = y/2/j,\ loge|, where /i is the rate of increase of the parameter and e is the 
amplitude of the noise, the field at last becomes 0(1) and the spatial pattern present 
is frozen in by the nonlinearity. Thereafter one observes spatial structure with char- 
acteristic size proportional to ( ' lo ^^ The increase of g need not be uniform, but 
the analytical estimates presented here rely on the rate of increase of g as it passes 
through being small. 
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Fig. 0.3. Mean number of zero crossings at g = 
and u = 0.01 ( uvve r line) and fi = 0.001 (lower line), 
simulations of (10. A). 



1. The solid lines are (0.21) with L = 2048 



The dots are obtained from finite difference 
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